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I. INTRODUCTION 

The need for a more fundamental theory at high en- 
ergies has been justified in many contexts. Divergences 
in quantum field theory, singularities in gravity and hav- 
ing a unified quantum framework for all forces, are some 
of them. A consequence arising from this consideration, 
which has been extensively studied, is the possibility of 
having Lorentz invariance violation given in the form of 
effective corrections [1]. This idea naturally leads to new 
extensions of the standard model and modified disper- 
sion relations for particles which may to be detected in 
sensitive experiments [2j. 

In this contexts the Myers-Pospelov theory is a model 
that introduces Lorentz invariance violation through di- 
mension five operators [H, Q ■ The breakdown of Lorentz 
symmetry that is included in the scalar, fermion and 
gauge sectors is characterized by an external timelike 
four-vector defining a preferred reference frame. Ex- 
perimental bounds for this model have been studied in 
several phenomena, such as synchrotron radiation Q , 
gamma ray bursts [f| , neutrino physics Q , radiative cor- 
rections Q . generic backgrounds [1, [l(| , and others [HI ■ 

In recent years and from a different approach, theories 
with higher time-derivatives have been proposed as an 
extension of the standard model of particles [l2[ • One of 
the main motivations is that these theories soften the ul- 
traviolet behavior of the quantum field theory, and hence 
problems like the hierarchy puzzle seems to be solved. 



Despite these theories contain negative norm states 
the theoretical consistency were established many years 
ago Although unitarity is maintained, the price to 
pay is the lost of microcausality [l5[ . 

The lost of microcausality and higher time derivatives 
are ingredients present in Myers-Pospelov theory. Then, 
it seems natural to wonder for the content of Lee- Wick 
modes in the Myers-Pospelov theory. In this work our 
main purpose is to rigorously establish these connection. 

The layout of this work is the following. In Sec. II 
we introduce the fermionic Myers and Pospelov model 
where we found the exact dispersion relation and we com- 
pute the corresponding eigenspinor in any background. 
In Sec. Ill we quantize the theory and we show explicitly 



the connection between the Myers-Pospelov model and 
Lee- Wick theories. In Sec. IV we discuss the origin of 
violations of microcausality and give the conclusions and 
final comments. 



II. FERMIONIC MYERS-POSPELOV MODEL 

The fermionic sector of the Myers-Pospelov theory is 
given by the Lagrangian 

£ = 4>(i@ - m) + 4(gi + 5275) (n ■ d) 2 yb, (1) 

where g\ and g 2 are inverse Planck mass dimension cou- 
plings constants and n is a dimensionless four vector 
defining a preferred reference frame with n 2 = +1,-1,0. 

Typically, phenomenological studies assume n purely 
in the temporal direction. However, in this section we 
will take n as general as possible and eventually we will 
consider some special choices. 

The variation of the Lagrangian ([T]) produces the equa- 
tions of motion 

[i$ - m + giiji{n ■ df + g2il b {n ■ d) 2 } ip(x) = 0. (2) 

In momentum space, 4>(x) — J d A pe~ w ' x ip(p), we obtain 
an algebraic equation, 



(</l - fa 5 - m) V> = 0, 



(3) 



13J where a M = - gin^(n ■ p) 2 and = g 2 n f _ l (n ■ p) 2 are 
four vectors which will help us to clear up the notation. 
Also, let us define the following matrices, 



£ = </l - fa 5 - 



(4) 



The first one, is related to the equation of motion. In fact, 
the solutions in the momentum space are those which sat- 
isfy Sip = 0, and then they are particular eigenvectors of 
£ . The second quantity is essentially the Pauli-Lubansky 
vector projected onto the fixed vector n, i.e., 



n-W = n^PpvMrp, 
= n^Pp v S ap , 



(5) 



2 



where M. and S are the total and spinorial Lorentz gen- 
erators, respectively. Indeed, 



2e^a p b u S, 



(J pi 

= 2g 2 (n-p) 2 e^p p n lJ S ap , 
= -2g 2 {n-p) 2 n-W. 



(6) 



To get this relation, we have made use of the fact that, 

[7",7"] 7 5 - l 2^ vap h«,l P ] = ^ vap S ap . 
These operators satisfy the interesting relations, 



£,h 



= 0, (£ + 2m)£ = a 2 - b 2 - m 2 - h. (7) 



This means that the solutions of the equation of motion 
can be expressed in terms of the eigenvectors of h, and 
that it must be satisfied the relation, 



m 



0, 



(8) 



which provides the general dispersion relation for the 
Myers- Pospelov theory, in terms of the eigenvalues h. By 
noticing that, 



h 2 = 4 [(a ■ b) 2 - a 2 b 2 ] I, 



(9) 



the eigenvalues of h are, h — ±2y/ (a ■ b) 2 — a 2 b 2 . Obvi- 
ously, this is true for (a ■ b) > ab. Otherwise, the operator 
h is antihermitian, and their eigenvalues are purely imag- 
inary with the same absolute value. Hence, for that case, 
we have a real dispersion relation given by, 

(p 2 -m 2 - 2 5l (n • pf + n 2 (g 2 - g 2 )(n ■ p) 4 ) * 

-4(n • pfg 2 2 ((n • pf - p 2 n 2 ) = 0. (10) 

Now, we will calculate the solutions of the equations of 
motion. As we pointed out above, we can find these 
solutions among the eigenvectors of h, hiph = hiph- Then, 
let us find those eigenvectors. To do so, we notice that 
the h operator can be written in terms of a rank two 
antisymmetric tensor, T^ v = a^b v — a^b^, that is, 



h = T^e^Svp. 



(11) 



From this tensor, we define two orthogonal three- vectors, 
(u) 1 = T° l = a°(bY -b°(aY =ue\, (12) 
(vT = \e ijk T jk - (a x bf = ve\, (13) 

and with them 

(w) % = (u x v) % = uve\, (14) 

where ei, e 2 and e% are three orthonormal space vectors 
on the direction of u, v and w, respectively. The norm of 
these vectors are, 

u = \/(a ) 2 (b) 2 + (b°) 2 (a) 2 - 2(a°&°)(o- 6), 



(a) 2 (6) 2 - (a-b) 2 . 



(15) 



Note that, 



T 2 = T^T^ = 2(v 2 -u 2 ), 

= 2(a 2 b 2 -(a-b) 2 ) = -h 2 



(16) 



The negative values of T 2 correspond to real eigenvalues 
for h and the positive ones correspond to pure imaginary 
eigenvalues. By making use of the analogy with the elec- 
tromagnetic tensor F we will call the T 2 < "electric" 
case and T 2 > the "magnetic" case. 

Now, we define the rotations and boosts generators in 
the spinor representation, 



(17) 



where the spatial indices are referred to the e basis de- 
fined above. Then, the h operator turns out to be 



h = -4(uJi + vJC 2 ) 



(18) 



Performing a non local boost-like transformation on the 
eigenspinor in the e 3 direction 



the h operator transforms as 



(19) 



h' = e ivK3 he- lr,K3 = -4 [(u cosh -q-vsmh rj)Ji 

+(ucoshry — usmh-q)^] . (20) 

Because — 1 < tanhr; < 1, we can distinguish two cases. 
For «>uwe can set tanh rt = - so that 



h' = -Ay/u 2 -v 2 J!. (21) 
However, for v > u, we can set tanh 77 = ^ such that, 

h' = -A^/v 2 -u 2 JC 2 . (22) 

Since the eigenvalues of J and K, are ± i and ± | respec- 
tively, we have h = ±2\/u 2 — v 2 = eh\h\ for u > v, and 
h = ±2Wv 2 — u 2 = ieh\h\ for v > u as we expected. 

The eigenspinors in the chiral representation for u > v 
can be written as 



(23) 



with (u- a)t;h = —ehu£h- Notice that these eigenvectors 
have the property, in the e basis, 

7Vfc = Cfc7°7Vfc- (24) 
However, for v > u, the eigenspinors have the form, 



IhXh 

Sh^Xh 



(25) 
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with (v ■ a)xh — £h v Xh- Similarly, the eigenspinors have 
the property, in the e basis, 



(26) 



The constants ah, fth, $h, lh reflects the fact that the 
eigenspinors are twofold degenerate. 

Now, we are ready to find the solutions of the equations 
of motion in terms of the spinors tph- Performing the 
same transformation on £ we obtain 

£' = e ir]K *£ e"^ 3 = $ - $' 7 5 - m, (27) 

where, a'^ = tv^a v and b'^ — A^6„ with A being some 
functions of the parameter 77. The non vanishing compo- 
nents of the four vectors a and b in the e basis are, 



a 1 

«3 



-{a- ei) = - 
-(a- e 3 ) = - 



a°(b-a)~b°(a)' 2 



(28) 



and, 



by = -(b-e 1 ) = — 



^{bf ~b°(a-b) 



b 3 = -(b-e 3 ) = -b -. 



(29) 



And the transformed components, in terms of the original 
ones are 

a' = a,Q cosh 77 + a 3 sinh 77, 

= ciq (cosh 77 sinh 77 

\ u 

a 3 = a 3 cosh 77 + ao sinh rj, 
= — ao y — cosh 77 — sinh 77 ) , 



(30) 



and, 



b' a = bo cosh 77 + 63 sinh 77, 

= &o ^cosh 77 sinh 77^ , 

63 = 63 cosh 77 + 60 sinh 77, 

= — 6q ( — cosh 77 — sinh 77 ) , 



(31) 



with the components along the directions e± and e 2 re- 
maining unchanged. 

In the electric case (u > v), we have a' 3 = b' 3 — and 



v 2 \h\ 



(32) 



In the other hand, in the magnetic case (v > u), we have 
a' = b' = and, 



iv 2 \h\ 
-ao\l -? - 1 = a — , 



(33) 



Hence, in the electric case the equations of motion, 
E'iph — 0, are 

[(a' - e h bx) 7 - {b' - Cfc oi) 7V - m] Vfc = 0, (34) 



where we have used (|24j) . This equation fixes the con- 
stants in the Eq. (|23|) 



Q<h = TVra, 

& = TV [(a' + 6^6!) - (b' + e h oi)] , (35) 

where TV is a normalization constant. The equations of 
motion, £'iph = in the magnetic case are 

[(01 - «e/ l 6 3 )7 1 - (61 - ie^a 3 )7 1 7 5 - m] i/j' h = 0, (36) 



where we have used property ()26[) . This implies that the 
constants in the Eq. ([2"5]) are 

7ft = TV' 771, 

S h = TV [(bi - ie h a' 3 ) - (a x - ie h b' 3 )] , (37) 

where TV' is another normalization constant. 

Finally, to get the complete solution to the general 
Myers-Pospelov equations of motion we must solve for 
p° = lu as a function of pi,P2,P3 in Eq. ((8|) or ([TO]) , 
substitute (gU), ([2"6). (|51|) and (j^J), and do the linear 
combination, 



t/j(x) 



7 . /; 



d 3 pa h (p)r h (p)e 



iojih (p)t—ip-x 



(38) 



where the index i runs over the solutions ujih(p) for a 
given ft, in the Eq. ([8]), and a l h (p) are integrations con- 
stants. To quantize the theory we must promote the con- 
stants a to creation and annihilation operators. Some 
peculiarities arising when no 7^ in the procedure of 
quantization will be discussed in the next section. In gen- 
eral, ((HJ is is a fourth order polynomial in uj for a given 
value of h, and it would yield at most four real solutions. 
The negative solutions correspond to antiparticles modes 
while the positive ones are particles modes. 

However, if n has no temporal component the degree 
of the polynomial in uj will turn to be of order two, and 
hence we will obtain a total of two solutions for each given 
value of h. This last situation is similar to the standard 
case where we have a total of four solutions. The former 
situation, namely, when no 7^ we are obtaining double 
number of solutions than in the standard case. This is 
due to the fact that we dealing with a theory with higher 
time derivatives as can be seen from the equation of mo- 
tion $2$ . In the next section we will discuss in more detail 
the nature of this extra solution. 

Let us Anally make two observations, the first one is 
that in the magnetic case we will not obtain real solutions 
for tu(p) because h has imaginary eigenvalues. The second 
one is that for n purely timelike h is proportional to the 
helicity and the energy will depend in general on this 
quantum number. Also we recover the well know results 
in this case [H, [l6| . 
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FIG. 1: The curves corresponding to the function f(u>) and 
the constant function g(uj) = p 2 . The intersection of both 
curves are the solutions of the dispersion relation (I39|l and 
Ei± are given in (|52|l . 



III. NON-LOCAL LEE- WICK MODES 

In the case n picks up a non vanishing temporal com- 
ponent the physics of the model shows up interesting fea- 
tures connected to the recently revisited model with in- 
definite metrics [lU, [l4| ■ For simplicity in this section 
let us assume a purely time vector n = (1,0,0,0) and 
cj2 = 0. The dispersion relation (|10l) reduces to 



f 



m 2 - 2gico 3 



(39) 



In order to see the qualitative behavior of the solutions let 
us define f(uj) = uo 2 — m 2 — 2giu> 3 + gfui 4 and g(uj) = p 2 , 
and plot these functions of lu in Fig. 1. 

The solutions are those points which intersects the 
curve / and the horizontal straight line in the upper 
half plane, which corresponds to the fixed value of the 
momentum square, i.e. g(u>) = p 2 . Hence, for small val- 
ues of p we find four solutions, one negative frequency 
which corresponds to an antiparticle and three positive 
frequencies. Among the positive frequencies the small- 
est one is the normal particle frequency and the other 
two correspond to Lee- Wick like modes, as we will show 
below. It is peculiar the behavior of one the Lee- Wick so- 
lutions whose frequency decreases with momentum, this 
will contin ue until th e momentum reaches the value of 
\p\ max ~ ~~ 771-2 wnere ^ collapses with the normal 

particle mode. Furthermore, it is worth noting the differ- 
ences in energy between particles and antiparticles which 
in the limit rag\ << 1 turns out to be 4 |<7i| m 2 . 

Let us discuss how the Lee- Wick modes came about. 
By looking at the equation of motion for the present case 
we see that this equation is a second order fermionic equa- 
tion. Then, let us consider a toy model which displays the 
essential properties. Let us take the (0 + l)-dimensional 
fcrmion model with higher time derivatives, 



The above Lagrangian can be rewritten as 

C = c±(id t — m+) c + — <y_{idt + m_) c_. (41) 
In terms of the new variables 

c + = a(iip + m^tp), c} + — a(—iip + ni-ip), (42) 

and, 

c_ = a(— iip — m + 'i[j), cL = a(iip — m + ?/>), (43) 

where a = { ? ] 2 an d 

Tl + VTTigm 



m± 



'Mi 



(44) 



Quantizing this model implies that the Hamiltonian is 

Hp = m + c^ + c + — m-C_C- — m_, (45) 

with the non vanishing anticommutation relations 

{c+,cl} = l, {c„,c t _} = -l. (46) 

This is a four state system with energies 0, m+, m_, m + + 
m_, and the eigenstates corresponding to the last two 
energies of negative norm states due to the wrong sign in 
the anticommutators above. 

It is worth noting that the standard first time deriva- 
tive is reached by performing the limit mg — > where 
m+ tends to m and m_ blows up to infinity as 1/g. The 
generalization to quantum field theory is straightforward 
and unitarity problems seems to be controlled by suitable 
definitions of propagators fl5| . 

In our case, recall the Lagrangian 

L = J d s x${i$ + 5 i7° dl - m)V>, (47) 
- J d 3 x ^ {id + 9l d 2 - h D )i), (48) 

where ho = id ■ V + m/3 is the standard Dirac Hamilto- 
nian operator. Now let us write the field in terms of the 
standard solutions of the Dirac Hamiltonian operator, 



s,i 



d 3 p 

TUT 



■.ut(p)^(p,t)e^ s , (49) 



Op 



L, 



-gip^ + iipip ~ mtpip. 



(40) 



where s is s spin index, i is the particle and antiparticle 
index, i.e., uf = u s and u| = v s in the standard notation 
and Eq p = \J p 2 + m 2 . Remembering that (p)it^ = 
2E 0p 5 sr 5 ij and hDuf^e^ = e.Eopuf { V )e^ s where 
ei = +1 and 62 = —1, we have 

L = Y J l d 3 pib?(P,t)(gid 2 o + id + e,£ 0p Mf (P,t).(50) 

s,i J 
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Now it is clear we have reduced the quantum field the- 
ory problem to a set of four quantum mechanical systems 
at a given momentum. Following our quantum mechani- 
cal example each of these systems can be decomposed in 
terms of two kinds of modes with indefinite metric, that 
is 

s,i 



where 



E l± = (Tl + y/l + tiAEkpgx) 



2<?i 



(52) 



It is worth noting that the fields i[>?± are non local in the 
sense that they are defined 



V4 = (td - E lT )t 



i±> 



(53) 



where Ei± are non local operators depending on the spa- 
tial derivatives. As in our quantum mechanical example 
the plus modes are standard particles with positive norm 
and the minus modes are Lee- Wick particles with nega- 
tive norm. Both modes have the degeneracy of spin and 
then we get a total of eight modes. There dispersion re- 
lation are those in (|52[) which agree with our qualitative 
argument in the beginning of the section. 



To conclude we would like to discuss microcausality in 
this model. Firstly, the fact that we are dealing with an 
explicitly Lorentz invariance violating model, we expect 
a small failure of microcausality at the kinematical level 
li [l8| . This can be checked by computing the first 
order corrections of the anticommutator of free fermionic 
fields at different spacetime points, 



{</>(*),#})} = / 



dp e 



Jc 



(a 2 (p) ~ b 2 (p) — m 2 ^ — h 2 (p) 



(54) 



where O is a local differential operator and C is a closed 
curve enclosing all the poles. However, there is another 
source of microcausality violation which come from the 
Lee- Wick nature of the model when the temporal compo- 

Even though this cannot been 



12] 



nent of n M is not zero 
seen from the anticommutator of the free fermions above 
we will obtain modifications after adding interactions to 
the theory. Then, the origin of the Lee- Wick failure of 
microcausality has a dynamical character. 

It would be interesting to study how this two sources 
of microcausality failure can interfere with each other. 
Finally we would like to mention that most of the results 
in this work could be generalized to the remaining sectors 
of the Myers-Pospelov theory. 



IV. DISCUSSIONS AND CONCLUSIONS 
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